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Abstract 

We  consider  a  cross-layer  optimization  problem  for  multi-hop  wireless  networks 
over  time  varying  channels.  The  system  consists  of  L  interfering  links,  where 
the  transmission  power  and  rate  of  link  l  (=  1,...,L)  are  specified  in  vectors 
P  =  [P\. . .  Pl\  and  X  =  [X\ . . .  Xl]  respectively.  In  every  time  slot,  the  scheduler 
schedules  the  transmissions  by  assigning  a  resource  allocation  vector  V  =  [P  X]T . 
We  denote  the  expectation  of  V  by  E(V).  Our  objective  is  to  find  the  optimal 
scheduling  policy  which  minimizes  the  cost  of  average  resource  consumption  while 
maintaining  average  service  guarantees  to  each  user.  We  develop  a  unified  frame¬ 
work,  in  which  the  cost  of  the  average  resource  consumption  is  given  by  a  convex 
function  f(E(V ))  and  the  minimum  average  service  guarantees  are  given  by  a  set 
of  convex  constraints  g(E(V))  <  0.  By  means  of  convex  optimization  and  sto¬ 
chastic  approximation,  we  obtain  the  solution  by  solving  the  corresponding  dual 
problem.  An  iterative  algorithm  is  proposed  and  analyzed,  which  schedules  the 
transmission  powers  and  rates  adapting  to  the  channel  variations.  If  the  channel 
states  are  described  by  a  finite-state  mixing  process,  it  is  shown  that  our  algorithm 
asymptotically  attains  the  optimal  cost. 


1  Introduction 

The  growing  interest  in  multi-hop  wireless  network  increases  the  demand  for  efficient 
resource  allocation  methods.  To  maintain  the  quality  of  service  and  manage  the  cost  of 
average  resource  consumption,  a  scheduler  is  responsible  for  coordinating  all  communi¬ 
cating  nodes  in  the  system.  In  this  paper  we  assume  that  the  scheduler  is  centralized. 
This  allows  us  to  obtain  insight  into  the  structure  of  the  problem,  and  our  solution  to 
the  problem  in  this  case  serves  as  a  benchmark  for  distributed  scheduling  algorithms. 

In  each  time  slot,  a  scheduler  specifies  the  transmission  power  and  rate  on  each 
link.  Since  wireless  channel  conditions  are  time  varying,  scheduling  a  transmission  over 
a  bad  channel  results  in  high  power  consumption  and  may  result  in  large  interference. 
Therefore,  an  intelligent  scheduler  should  take  into  account  the  instantaneous  channel 
state  information  (CSI),  such  as  background  noise,  pathloss  gain,  and  fading,  as  well  the 
variations  of  these  quantities  in  time  and  space,  in  order  to  achieve  an  efficient  resource 
utilization. 


Several  studies  have  addressed  the  problem  of  throughput  utility  maximization,  con¬ 
sidering  time  varying  channels.  In  [9],  Tse  proposed  the  proportional  fair  scheduler  (PFS), 
which  exploits  multi-user  diversity  in  a  time  varying  environment.  Under  the  stationary 
assumptions  on  the  channels,  the  throughput  vector  [XfV9, . . . ,  Xff9]  of  PFS  achieves  the 
maximum  of  the  logarithmic  utility  function  y);=1  log(A"“1,s)  among  all  schedulers  [5]. 
Agrawal  [1]  and  Stolyer  [8]  further  generalize  this  idea  by  considering  arbitrary  utility 
functions  of  the  form  i  /i(A“1'9),  where  /;(•)  is  concave  and  differentiable. 

One  weakness  of  the  utility  maximizing  scheduling  is  that  the  optimal  scheduler  favors 
users  with  better  channel  conditions.  This  implicitly  decreases  the  priority  of  users  with 
poor  channel  quality.  However,  it  may  be  highly  desirable  to  impose  a  minimum  average 
performance  requirement  on  each  link.  To  address  this  issue,  Lee  et  al.  [7]  reconsidered 
the  throughput  maximization  problem  and  impose  a  minimum  throughput  constraint 
on  each  user.  However,  their  solution  is  only  asymptotically  optimal  in  the  limit  of  an 
infinite  number  of  users. 

The  structure  of  optimal  scheduler  under  time  varying  channels  is  investigated  in  this 
paper,  in  which  we  minimize  the  cost  of  average  resource  consumption  while  maintaining 
a  minimum  average  performance  guarantee  on  each  link.  The  problem  is  formulated  as 
an  optimization.  An  iterative  algorithm  is  proposed  to  obtain  the  solution.  By  means 
of  convex  optimization  and  stochastic  approximation,  we  verify  that  our  algorithm  ap¬ 
proaches  the  optimum  irrespective  of  the  number  of  users.  We  show  that  besides  the  link 
layer  scheduling  and  power  control  problem,  our  framework  can  be  further  extended  to 
solve  other  network  layer  optimization  problems,  such  as  power  efficient  routing. 

2  Notations  and  System  Model 

We  consider  an  omni-antenna  wireless  system  consisting  of  N  nodes  and  L  links.  The 
instantaneous  transmission  power  and  rate  are  specified  in  the  vectors  P  =  [P\ . . .  Pl \ 
and  X  =  [Ai  . . .  Xj]  respectively.  Time  is  divided  into  fixed  duration  intervals,  called 
time  slots.  At  the  beginning  of  time  slot  k  (=  l,...,oo),  the  scheduler  specifies  the 
resource  allocation  for  the  system  in  the  vector  V  =  [P  X]T  .  We  model  the  underlying 
time  varying  phenomena  as  a  stationary  process  {co(k),k  >  0}  with  finite  state  space 
S  and  stationary  distribution  {s(cu),ct;  G  5}.  The  outcome  of  the  process  is  called  the 
system  state.  In  state  ui,  the  power  loss  gain  from  the  transmitter  of  link  l\  to  the 
receiver  of  link  l2  is  denoted  by  including  path  loss  and  fading.  On  link  /,  the 

receiving  end  experiences  noise  power  ffco).  We  hence  have  the  signal  to  interference 
and  noise  ratio  (SINR)  7 fu)  =  ^ •  The  sets  of  all  outgoing  links  and 

incoming  links  at  node  n  are  denoted  by  £(n)  and  P{n)  respectively.  Assuming  the 
maximum  instantaneous  data  rate  possible  by  link  l  is  a  function  R^i)  of  SINR  and  the 
peak  transmission  power  of  node  n  is  limited  to  Pffax ,  a  resource  allocation  V  is  feasible 
if  and  only  if  X  and  P  satisfy  the  following  constraints. 

(Power  Constraints) 

Pi  >  0  for  l  —  1, . . . ,  L  (1) 

Elesn  pi  <  Pnax  for  n  =  1, . ,  N 


(Rate  Constraints) 


A;  >  0  for  l  —  1, . . . , L 

A,  <i?(7,H)  for  1  =  1,...,  L 


(2) 


We  denote  T>{u)  as  the  set  of  all  feasible  resource  allocation  V.  It  is  uniformly  bounded 
over  S  under  the  finite  state  assumption. 


3  Scheduling  Policy  and  Achievable  Set  of  Average 
Resource  Allocation 

In  this  paper,  we  consider  a  probabilistic  scheduling  policy,  where  the  scheduling  deci¬ 
sions  are  random  variables  over  the  state  space  S.  Specifically,  in  state  oj,  a  scheduling 
policy  it  specifies  a  feasible  resource  allocation  V  randomly  from  a  specified  control  set 
JC(u)  with  probability  distribution  by,(V),  where  JC(u)  is  a  subset  of  T>(u),  (f>w(V)  >  0 
and  ’Yhv&K(u)  A(^)  —  1-  Therefore,  given  a  scheduling  policy  7 r,  the  average  resource 
allocation  En(V)  can  be  evaluated  as  below. 

E,(V)  =  Y.  UV)V.  (3) 

VgK(uj) 

We  omit  the  subscript  n  when  the  policy  is  known  explicitly.  We  define  the  achievable 
set  V  =  {hj  There  exists  it  such  that  v  =  En(V) j  .  The  uniform  boundedness  of  T>{u) 

implies  the  boundedness  of  D.  Moreover,  one  can  verify  the  convexity  of  T>  by  a  time¬ 
sharing  argument.  In  the  next  section,  we  introduce  the  framework  for  optimal  link 
scheduling  and  power  control,  where  we  only  look  into  link  layer  scheduling  problems. 
The  extensions  to  network  layer  routing  problems  are  discussed  in  section  5. 

4  Generalized  Framework  for  Link  Scheduling  and 
Power  Control 

Considering  link  scheduling  and  power  control,  most  of  the  link  layer  resource  manage¬ 
ment  problems  can  be  formulated  as  a  special  case  of  the  following  optimization, 

minimize  f(E(V))  (4) 

subject  to  g(E(V))  <  0 
E(V)  e  v, 

where  /(•)  is  a  scalar  valued  convex  function  and  g(-)  is  a  vector  valued  convex  function. 
For  example,  if  we  set  the  cost  function  to  f(E(V))  =  —  ]C;=|  log(i?(W)),  the  optimiza¬ 
tion  reduces  to  the  proportional  fair  scheduling  problem.  If  we  let  f(E(V))  =  i  -^(P/), 
it  becomes  the  power  efficient  scheduling  problem. 

Since  the  constraint  set  T>  is  bounded  and  convex,  if  E(V)  is  considered  as  the  in¬ 
dependent  variable,  the  problem  (4)  fits  in  with  the  framework  of  convex  optimization. 
The  duality  technique  [2]  in  convex  analysis  helps  to  simplify  the  problem  by  moving  the 
constraints  into  the  cost  function.  The  dual  function  q(/3)  is  defined  as 

q0)  =  min  L(E(V)J),  (5) 

E(v)e  v 

where  L(E(V),  j3)  :=  f(E(V))  +  (3  ■  g(E(V))  is  the  Lagrangian  of  the  primal  problem. 
The  dot  product  is  defined  as  the  inner  product  of  two  vectors.  One  can  verify  that  q(/3) 


is  a  concave  function  [2].  The  duality  theory  defines  the  dual  problem  below. 

maximize  q(P)  (6) 

subject  to  P  >  0 

The  solution  of  (6)  is  denoted  by  P*,  where  the  elements  of  the  vector  /3*  are  called 
Lagrange  multipliers.  From  the  weak  duality  theorem,  it  is  recognized  that  the  dual 
optimal  value  q*  is  a  lower  bound  of  the  primal  optimal  value  f*.  If  there  exists  a  feasible 
solution  in  the  relative  interior  of  T> ,  the  strong  duality  theorem  ensures  that  q*  =  f* . 
In  the  following  subsections  we  first  solve  the  problem  with  a  linear  cost  function  and 
linear  constraints,  and  then  generalize  it  to  nonlinear  problems  with  the  help  of  additional 
auxiliary  variables. 


4.1  Linear  Cost  and  Linear  Constraints 


For  a  system  with  a  linear  cost  function  f(v)  =  a  ■  v  and  constraints  g(v )  =  An  —  b , 
the  primal  problem  (4)  becomes 

minimize  a  ■  E(V)  (7) 

subject  to  A E(V)  <  b 

E(V)  e  v, 

where  a  and  b  are  vectors  and  A  is  a  matrix. 

Theorem  1.  Given  the  optimal  dual  solution  P* ,  the  optimal  scheduling  policy  i r*  for  (7) 
is  given  by  {  ( 0*  (V)  )  |  u>  G  S  },  the  optimal  control  sets  and  the  corresponding 
distributions,  where  JC*(u)  =  argmin^eX)^  a  ■  V  +  P*  ■  ( AV  —  b)  ,  and  0*(V’)  is  the 
distribution  over  such  that  E7T( V)  e  V,  g(EK(V))  <  0,  and  /3*  ■  g(En(V))  =  0. 

Proof.  The  dual  function  of  (7)  can  be  written  as 


Q  (P)  = 


=  mm 

E(V)£T> 


a  ■  E{V)  +p-(AE(V) 


(8) 


(Oi 


=  nun  >  stcn)  ^  PuiY)  «  •  V  +  P  ■  (AV  —  b) 

vgK(uj) 

(c 0)  min  {  Pcj(V)  a  ■  V  V  P  ■  (AV -b)  } 

K{u),<t>„(V)  -  L  J 

^e*s  V£K.{uj) 


To  achieve  the  minimum  in  (9),  the  control  set  JC(u)  must  be  a  subset  of 


arg  mm 
vev(u>) 


a-V  +  p-  (AV  -  b) 


(9) 

(10) 

J 

The  the  proof  is  accomplished  by  the  following  lemma.  □ 

Lemma  1.  (Optimality  Conditions  for  Convex  Optimization)  [2].  The  primal  and  dual 
variable  pair,  (E(V*),  P*),  is  an  optimal  solution-Lagrange  multiplier  pair  of  (f)  if  and 
only  if 

E(V*)  G  T>,  g(E(V*))  <  0  (Primal  Feasibility) 

P*  >  0,  (Dual  Feasibility) 

E{y*)  —  argminE^e:p  L(E(V),  P*)  (Lagrangian  Optimality) 

P*  ■  g(E(V *))  =  0  (Complementary  Slackness) 


Since  {4>u (V)}  sum  up  to  one,  from  (9),  we  can  further  deduce  that 


\  s(u>)  min 


E 


a  ■  V  +  0  ■  (AV  -  b) 


min  a  ■  V  +  (3  ■  ( AV 
lvev(u) 


For  ease  of  exposition,  we  define  the  quasi-dual  function  as 


q{uj) 


min  <  a  ■  V  +  (3  ■  ( AV  - 

vev(u>)  f 

a-V*(u,0)+P-  (AV*  (u,  P)-b), 


(ii) 


(12) 

(13) 


where  V*(u,/3)  is  a  solution  to  the  minimization  on  the  right  hand  side  of  (12).  Simple 
convexity  arguments  ensure  that  q(co,/3)  is  concave  in  (3.  Substituting  q(u,/3)  into  (11) 
and  (6),  the  dual  problem  is  turned  into  the  following  stochastic  convex  optimization 
problem. 

maximize  E[q(u,/3)\  (14) 

subject  to  (3  >  0 

For  any  concave  function  q(/3),  a  vector  jl  is  called  a  subgradient  of  q(/3)  if  it  satisfies 
the  inequality  q(/3')  —  q(j3)  <  /2  -  (/?'  —  /?)  for  all  (3 By  Danskin’s  theorem  [2],  one  can 
show  that  the  vector  (A V*(u,(3)  —  b)  is  a  subgradient  of  q(u,(3)  with  respect  to  (3 . 

To  solve  (14),  we  apply  the  following  algorithm. 

Algorithm  1.  (Quasi- Gradient  Method  [3])  At  the  beginning  of  the  kth  time  slot,  update 
the  dual  variables  (3k  by  the  recursion  (15)  below  with  step  size  €k,  where  ek  >  0,  )T)fc  = 
og,  and  Yjk  4  <  °° . 


P  =  [Pk~l  +  ^k-ijdk~l}+ ,  (f3°  =  0,  jf  —  0)  (15) 

The  vector  jlk  is  a  subgradient  of  q(u(k),Pk)  and  [  ■  ]+  :=  max(-,0).  We  choose  the 
vector  jlk  =  ( AV*(k )  —  b)  in  particular,  where 


V*(k)  G  JC*(u>(k)) 


arg  mm 
vev(u>(k)) 


a 


V  +  (3k  ■  (AV 


If  {u(k),k  >  0}  are  independent  identically  distributed  (i.i.d.),  the  algorithm  ensures 
that  (3k  — >  (3*  as  k  — >  oo  with  probability  one.  In  addition,  if  we  choose  the  step  size 
efc  =  o>k/k,  where  linq,  a/..  — >  a  and  a  >  0,  the  long-term  average  of  V*(k)  converges  to 
the  primal  optimal  solution  of  (7)  with  probability  one. 

Remark  1.  The  i.i.d.  assumption  on  u(k)  can  be  weaker.  Algorithm  1  is  still  applica¬ 
ble  if  the  dependency  among  {u(k),k  >  0}  decreases  in  time.  When  the  dual  solu¬ 
tion  (3*  is  unique,  one  can  verify  via  the  stochastic  approximation  (  [4],  [6]  Theorem 
8.2.5  and  5.2.2)  that  (3k  — ■>  (3*  with  probability  one  if  u(k)  satisfies  the  mixing  property: 
\Ekpt,((3,oj(i))  —  Efi((3,u>(i))\  — >  0  as  i  ^  oo,  where  Ek  denotes  the  expectation  condi¬ 
tioned  on  Ek,  the  filtration  of  {/31,  ca(l), . . . ,  j3k ,u>(k)}. 


To  prove  the  asymptotic  optimality,  we  rewrite  (15)  as 


=  Pk~l  +  j +  i*"1)  ,  (16) 

where  5*  is  the  vector  of  reflection  terms  of  the  [  •  ]+  operator  and  i*  >  0.  Considering 
a  telescoping  sum  on  both  sides,  we  derive  the  following. 

Pf  =  Pt-1+^(atpt-1+atzt-1+pt~1-pt-1)  (17) 

1  t-i  t-  i  t-i  t-i  t-i 

=  ^  -  a^k  +  f^2^ak~ +  ~t  a$k  +  7  +  7  S 

/c=0  fc=0  /c=0  /c=0  k= 0 


Since  (3k  — >  /?*  with  probability  one,  we  have  lim^oo  j  ^jj.=Q  Pk  =  /3*  with  probability 
one.  The  boundedness  of  V(u)  implies  that  jlk  =  (A V*(k)  —  b )  is  bounded.  Because  7* 
is  the  reflection  term,  it  is  also  bounded.  Note  that  lim k^ooi^k  ~  a)  =  0.  Taking  limits 
on  (17)  and  then  subtracting  /3*  from  both  sides,  we  arrive  at 


lim 

k — »oo 


0  with  prob.  1 


(18) 


When  (3^,  the  mth  entry  of  P*,  is  positive,  the  corresponding  reflection  term  z(n 
vanishes  for  k  large  enough.  By  (18),  we  therefore  have  lini/,.^oc  |  Y2t=o  fAn  =  0-  Similarly, 
if  fi*n  =  0,  one  can  show  that  lim  sup^^  \  Ylt= o  dm  <  0.  Replacing  jlk  with  (A V*(k)—b), 
we  get  the  asymptotic  feasibility  described  in  the  theorem  below. 

Theorem  2.  (Feasibility)  If  the  limiting  time  average  of  V*(k)  exists,  it  lies  in  the 
feasible  region  of  (7)  with  probability  one,  which  means 

—  b  <  0  with  prob.  1. 


1  *_1 

lim  Tr(fc) 

t — XX)  t  ‘  ^ 

k= 0 


To  show  the  asymptotic  optimality,  we  take  time  average  on  (13)  and  then  group 
the  summands  with  respect  to  system  state  a The  theorem  below  is  then  derived  by 
applying  the  stationary  property  of  ou(k). 

Theorem  3.  (Optimality)  If  the  limiting  time  average  of  V*(k)  exists,  it  attains  the 
optimal  solution  with  probability  one.  That  is 


r 


with  prob.  1. 


4.2  Convex  Cost  or  Convex  Constraints 

If  /(•)  or  g(-)  are  nonlinear  functions,  we  cannot  apply  the  decomposition  in  (9)  directly 
using  the  separable  property.  To  accommodate  the  non-separable  issue,  we  define  an 


auxiliary  variable  Y  to  substitute  E(V)  in  /  and  g.  Assuming  the  elements  of  T>  are 
bounded  above  componentwise  by  a  vector  Drnax ,  the  optimization  (4)  becomes 

minimize  /(F)  (19) 

subject  to  g(Y)  <0,  0  <  Y  <  Dmax 

Y  =  E(V),  E(V)  G  V 

By  weak  duality  theory  and  convexity,  the  solution  can  be  approached  using  the  algorithm 
similar  to  that  for  the  linear  cases.  Note  that  if  Dmax  is  not  known  explicitly,  one  can 
use  an  estimate  instead  and  update  it  in  each  time  slot. 


5  Optimal  Routing  and  Flow  Rate  Assignments 

Now  we  extend  our  framework  to  consider  other  network  layer  optimizations.  There  are 
J  flows  in  the  system  indexed  by  j  —  1, . . . ,  J.  The  source  and  destination  nodes  of  flow 
j  are  denoted  by  ns(j)  and  nd(j)  respectively.  All  links  used  by  flow  j  for  routing  form 
a  set  G(j)  of  size  rrij.  The  collection  of  all  flows  on  link  l  constitutes  the  set  7i(l).  In  a 
time  slot,  link  l  has  data  rate  C\j  for  flow  j  G  lH.il).  Given  G(j)  =  {h,  -  •  - ,  we  define 
Cj  =  [Cia  . . .  Ci  j]  as  the  routing  vector  of  flow  j.  Furthermore,  we  group  all  routing 

vectors  into  a  vector,  C  —  [Cj  ...  Cj\ ,  called  the  flow  vector.  The  essential  requirement 
for  routing  problem  is  to  maintain  minimum  end-to-end  throughput  Cj  on  each  flow  j. 
Under  this  requirement,  two  constraints  are  imposed  on  { Cij }. 


(Flow  Constraints) 


^jeHii)  (  'J  —  X  ’  ^  —  • 

ctj  >  0 


,,L 


(20) 


(Flow  Conservation  Constraints) 

For  n  —  1, . . . ,  N.  j  =  1, . . . ,  J 

E  E(Clj)  -  E  +  Vn3  =  wllere  vnj 

ie£(n),  ieg(j)  i&F(n),  ieg(j) 


Cj,  if  n  =  nd(j) 

-Cj,  if  n  =  ns(j)  (21) 
0,  otherwise 


The  resource  allocation  vector  now  becomes  V  =  [?  I  C]T.  In  particular,  the  new  feasi¬ 
ble  set  consists  of  all  vectors  V  —  [P  X  C]T  satisfying  the  power  constraints,  rate 

constraints  and  flow  constraints.  In  the  following  subsection,  we  discuss  our  framework 
for  the  power  efficient  routing  problem. 


5.1  Power  Efficient  Routing 

The  power  efficient  routing  problems  can  be  formulated  as  the  optimization  below. 

L 

5>W)  (22) 

1=1 

E[V]  satisfies  (21) 

E(V)  G  V 


minimize 


subject  to 


The  quasi-dual  function  of  this  problem  is  given  below. 

L  J  N 

nun  +  Cli~  °ij  +  unj)  (23) 

Ve T>(u)  i=l  j=\  n=\  le£(n),leG(j)  l&T(n),  ieG(j) 

The  dual  problem  is  therefore  expressed  as 

maximize  Eq(u,/3) 

subject  to  0  E  Mm,  (24) 

where  m  is  the  dimension  of  the  constraints.  To  evaluate  (23),  first,  we  rearrange  the 
summands  on  the  right  hand  side  of  (23)  with  respect  to  Cy, 

L  J  L  J  N 

q(u,/3)  =  _min  ^  Pi  +  EE  aijCij  +  (T0,  where  cr0  =  EE  fikj  V nj  (25) 

V&V{u)  l=1  j=1  l=1  j=i  n=  1 

and  uij  is  the  new  coefficient  of  Cij  after  rearrangement.  For  fixed  P  and  X,  equation 
(25)  is  a  linear  optimization  problem  in  Cij  and  the  flow  constraints  form  a  polytope  for 
{Cij}.  Hence  it  has  solutions  at  extreme  points  of  (20).  One  can  verify  that  the  optimal 
flow  rates  borne  on  each  link  are  given  by 

{Xi  if  aij  =  argminieW(p  atj,  and  atj  <  0 

(If  not  unique,  pick  up  one  randomly)  (26) 

0  else. 


Substituting  (26)  into  (25)  and  rearranging  the  summands  according  to  Xi,  we  arrive  at 
the  equation 

L  L 

q(ui,  0)  =  min  "jP  Pi  +  ^  +  a0,  (27) 

v^o)tr  tr 

where  a[  is  the  corresponding  coefficient  after  the  second  rearrangement.  Similarly,  fixing 
P,  the  optimal  transmission  rate  on  link  /  is  given  by 

Ri(li)  if  cr[  <  0 
0  if  a[  >  0. 

Substituting  (28)  into  (27)  and  then  grouping  the  summands  according  to  Ri^i),  the 
evaluation  of  quasi-dual  function  reduces  to  the  following  optimization  with  respect  to 
transmission  powers: 

L  L 

q{v,  P)  =  min  TR  +  T  cr'/RM  +  0o, 

where  a'/  is  the  corresponding  coefficient.  The  solution  of  (22)  then  can  be  obtained 
through  the  following  algorithm. 

Algorithm  2.  At  time  slot  k,  the  optimal  scheduling  policy  schedules  a  transmission 
power  vector  P*(k)  from  the  set  below. 

L  L 

E^+E^'W 

i=i  i=i 


arg  mm 

V£V{uj{k)) 


(29) 


The  optimal  data  rate  X*(k)  and  flow  rate  Cf(k)  can  be  deduced  from  (28)  and  (26). 
The  dual  variable  is  updated  recursively  as 

J2  CZ (k)  -  J2  G’i  W  +  "nj  ,  (  P“  =  o  ) 

l££(n)  l£^F(n) 

where  lim^oo  =  a  and  a  >  0. 

Let  V*(k)  =  [ P*(k )  X*(k )  C*(k)}T,  if  the  process  Lv(k)  is  i.i.d.  or  mixing  as  described 
before,  the  long-term  average  of  V*{k )  converges  to  the  optimal  solution  of  (22). 

In  the  next  section,  we  exam  a  numerical  example  of  power  efficient  routing,  where 
we  apply  the  linear  rate  function  R{pi)  =  W'ji.  The  cost  function  in  the  evaluation  of 
(29)  then  becomes  componentwise  concave  in  P One  can  show  that  the  solutions  to 
such  problems  happen  at  the  extreme  points  of  power  constraints  (1). 

6  Numerical  Examples 

This  example  consists  of  7  nodes  and  8  links,  which  are  depicted  in  figure  1(a)  .  There 
are  two  flows,  each  requires  minimum  throughput  C.  Flow  1  originates  at  node  1  and  is 
destined  to  node  5;  it  exploits  links  1,  3,  4,  5,  and  7  to  route  the  traffic.  Flow  2  originates 
at  node  3  and  is  destined  to  node  7;  it  uses  link  2,  3,  4,  6,  and  8  to  route  the  traffic.  We  set 
pmax  cjQ  rnyy  ancj  \y  50  MHz.  The  channel  states  are  i.i.d.  The  background  noises 
are  normally  distributed  with  mean  zero  and  an  average  power  of  1  mW.  Moreover,  the 
noise  power  is  truncated  above  at  10  mW  and  rounded  to  the  second  digit  in  precision. 
The  channel  gain  G/p2  is  given  by  e/d2  (A,  12),  where  e  is  an  exponential  random  variable 
with  mean  1  and  d(l\,  )  is  the  distance  between  the  transmitter  of  link  li  and  the  receiver 

of  link  Z2-  In  addition,  the  factor  e  is  truncated  from  above  at  2  and  rounded  to  the  third 
digit  in  precision.  The  step  size  is  =  5020(/fc .  To  investigate  how  system  performs 
with  respect  to  C,  we  gradually  increase  the  minimum  throughput  requirements  from 
0.5  to  30  Mbps.  The  traces  of  the  average  flow  rates  carried  on  each  link  are  plotted 
in  figure  1(c).  The  upper  curve  contains  8  overlapped  traces  of  the  time  averages  of 
{C11,  C22,  C31,  C32,  C41,  C42,  C51,  C62}  and  the  lower  curve  2  overlapped  traces  of  the  time 
averages  of  {C71,  C$2} ■  Because  the  topology  is  symmetric,  we  only  focus  on  the  behavior 
of  flow  1.  It  is  recognized  that  when  the  requested  throughput  is  below  5  Mbps,  the 
optimal  route  of  flow  1  is  { 1  — >2— >4^6— >5};  as  the  requested  throughput  increases, 
the  direct  route  {1  — >  5}  comes  into  play.  Although  both  {1  — ^  2  — ^  4  — ^  6  — ^  5}  and 
{1  — »  5}  are  energy  efficient  paths,  on  route  {1  — >  2  — ^  4  — ^  6  — ^  5}  each  link  contributes 
less  interference  to  the  system,  and  the  possibility  that  all  of  the  links  are  in  deep  fade 
at  the  same  time  is  small.  In  figure  1(b),  we  gather  the  data  from  slot  5001  to  slot  25000 
and  plot  the  average  power  consumption  and  average  quasi-dual  value  respect  to  C .  We 
can  see  that  the  curve  of  average  power  consumption  closely  follows  that  of  the  average 
fl]f)  as  expected.  The  gap  between  both  curves  will  diminish  to  zero  in  the  limit  of 
an  infinite  number  of  time  slots.  When  C  is  small,  the  scheduler  works  like  a  TDMA 
system.  In  other  words,  in  every  slot,  at  most  one  of  the  links  is  activated,  and  the 
scheduler  refrains  from  transmitting  on  a  link  unless  that  channel  is  in  good  condition. 
Therefore,  the  total  power  consumption  increases  linearly  in  the  low  throughput  region. 
However,  as  C  goes  up  to  20Mbps,  more  links  need  to  be  activated  in  a  slot.  The  total 
power  consumption  then  increases  nonlinearly  due  to  interference. 


(a)  Network  Topology 


(b)  Power  Consumption 


(c)  Average  Flow  Data  Rates 


Figure  1: 


7  Conclusions 


In  this  paper,  we  have  shown  how  a  joint  link  scheduling,  power  control,  and  routing 
optimization  problem  over  time  varying  channels  can  be  formulated  as  a  stochastic  convex 
optimization.  In  addition,  we  proposed  a  centralized  iterative  algorithm  which  exploits 
the  CSI  to  schedules  the  transmissions.  We  have  proved  that  the  optimality  and  feasibility 
is  attained  asymptotically.  The  numerical  results  show  that  for  flows  with  minimum 
throughput  requirements,  to  save  power  consumption  the  scheduler  opportunistically 
routes  the  traffic  through  links  with  better  channel  conditions. 
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